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ABSTRACT 

We report the validation and characterization of three new transiting exoplanets using SOPHIE radial velocities: KOI-614b, KOI-206b, 
and KOI-680b. KOI-614b has a mass of 2.86 + 0.35 M Jup and a radius of 1.13^ 2 g R Jup , and it orbits a GO, metallic ([Fe/H]=0.35+0.15) 
dwarf in 12.9 days. Its mass and radius are familiar and compatible with standard planetary evolution models, so it is one of the few 
known transiting planets in this mass range to have an orbital period over ten days. With an equilibrium temperature of T eq - 
1000 + 45 K, this places KOI-614b at the transition between what is usually referred to as “hot” and “warm” Jupiters. KOI-206b 
has a mass of 2.82 + 0.52 M Jup and a radius of 1.45 + 0.16 R Jup , and it orbits a slightly evolved F7-type star in a 5.3-day orbit. It 
is a massive inflated hot Jupiter that is particularly challenging for planetary models because it requires unusually large amounts of 
additional dissipated energy in the planet. On the other hand, KOI-680b has a much lower mass of 0.84 + 0.15 M Jup and requires less 
extra-dissipation to explain its uncommonly large radius of 1.99 + 0.18 R Jup . It is one of the biggest transiting planets characterized so 
far, and it orbits a subgiant F9-star well on its way to the red giant stage, with an orbital period of 8.6 days. With host stars of masses 
of 1.46 + 0.17 M q and 1.54 + 0.09 M 0 , respectively, KOI-206b, and KOI-680b are interesting objects for theories of formation and 
survival of short-period planets around stars more massive than the Sun. For those two targets, we also find signs of a possible distant 
additional companion in the system. 
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1. Introduction 


Since 2010, we have been conducting a radial-velocity follow¬ 
up of Kepler (Borucki et al. 2010) tra nsiting planetary candi¬ 
dates with the SO PHIE spectrograph ( [Perruchot et al.|[2008[ 
|Bouchy et al.|2013| ). Transits combined with radial velocity mea¬ 
surements enable us to infer the masses and radii of the star 
and p lanet up to a one-parameter degeneracy ( Wright & Gaudi 
|2013| ). To derive the absolute planetary parameters, an indepen¬ 
dent stellar parameter is needed, which can be estimated with the 
aid of stellar models. With the mass and radius of the planet, we 
can determine the mean density that gives insight into the com¬ 
position and formation of the objects. In addition, such a survey 
allows the rate of false-positive of the Kepler close-in giant exo¬ 
planet candidates to measured (Santeme et al. 2012). 

In this paper, we present the validation and characterization 
of three new transiting exoplanets: KOI-614b, KOI-206b, and 


* Based on observations made with SOPHIE on the 1.93-m telescope 
at the Observatoire de Haute-Provence (CNRS), France. 


KOI-680b. Two of them, KOI-206b and KOI-680b, have already 
been announced by o ur group to be planetary companions in 
|Santerne et al.| ( [20T2| ). The third one, KOI-614b, was not in the 
sample considered by |Santerne et al. ( |2012| ) owing to its transit 
depth of 0.39%, slightly lower than the 0.4% selection cut-off. 

The structure of the paper is as follows. We present the 
Kepler and SOPHIE observations (Section [2]), the derivation of 
the spectral parameters (Section [3]), the stellar and planetary pa¬ 
rameter estimation (Section [4]), the planetary evolution models 
(Section [5}, and finally the discussion (Section [6]) and the con¬ 
clusions (Section [ 7 ]) of our work. 


2. Observations 

2.1. Photometric detection with Kepler 

The three targets analyzed in this paper were identified as KOI 
(Kepler Object of Interest) by |Borucki et al.|p011| ). Their IDs, 
coordinates, and magnitudes are reported in Table [2] Two of 
them were found to be hosting a transiting companion of plane- 
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tary nature with high probability, namely KOI-206 and KOI-680 
(vetting flag of 2), whereas KOI-614 was considered a moder¬ 
ate probability candidate because of the V-shaped transit (vet¬ 
ting flag of 3). There is no significant spatial shift at the time 
of the transits in the Kepler photometric centroid for the three 
targets. No other transits with different periods were detected in 
any of the light curves, so there are no signs of multiple tran¬ 
siting systems. At the time of the beginning of the analysis, the 
Kepler light curves of the first 15 quarters are publicly available 
at the MAST archiv^] Long cadence (29.4 minutes sampling) 
data only is available for KOI-614, but KOI-206 and KOI-680 
have short-cadence data (1-minute sampling) available in 3 of 
the 15 first quarters. We used the light curve of quarters Q1 
to Q15 calibrated through the Kepler Science Pipeline and cor¬ 
rected for instrumental artifacts and other known systematics er¬ 
rors ( [Jenkins et al.|2010) . The three light curves present transits 
with depths of about 0.5% and a typical uncertainty for individ¬ 
ual points at the level of 160-280 ppm for the long-cadence data. 

Besides the transits, the three light curves show signs of 
stellar variability with flux modulations of about 90, 75, and 
45 ppm amplitudes for KOI-206, KOI-614, and KOI-680, re¬ 
spectively. As a_result, we computed their Lomb-Scargle peri- 
odograms (Press & Rybicki|1989| using the long-cadence sam¬ 
pling after subtracting the transits (Fig. |T] top panels). They all 
present peaks compatible with the rotational period of the stars 
as expected from their measured projected rotational velocities 
(v sin /*, Section]?} and radii (Section]?]), which give upper limits 
at 10.5 ± 1.5, 23.3L£, and 27.3 ± 5.7 days for KOI-206, KOI- 
614, and KOI-680, respectively. Given the amplitude of the vari¬ 
ations and their period, we conclude that those peaks are likely 
due to spots covering the stellar surface. The evolution of the 
spots, which vary in number and position with time, may scram¬ 
ble the periodicity of the signal, so we divided the light curves 
in chunks of the same duration, chosen to span over about five 
expected rotation period, and we again computed the Lomb- 
Scargle in each one. For each light curve, the resulting mean 
periodograms were plotted in Fig. [l] (middle panels). We also 
give the autocorrelation functions computed on the whole light 
curve (Fig. |T] bottom panels). 


The main period and its harmonics are clearly detectable for 
the three targets. All methods are in good agreement, and we 
estimate stellar rotational periods of 7.95 ± 0.90, 18.09 ± 0.10, 
and 28.24 ± 0.23 days for KOI-206, KOI-614, and KOI-680, 
respectively. Finally, owing to the high level of noise at high 
frequency, we cannot detect any evidence of p ulsations in the 
expected range for those stars ( Belkacem et al.|20TT) . 


2.2. Radial velocities with SOPHIE 

The spectroscopic follow-up of the three targets was performed 
with the SOPHIE spectrograph (Perruchot et al. 2008| ) after its 
upgrade made in 2011 ( [Bouchy et al.||2013| . The three Kepler 
targets were observed from July 2011 to August 2012 with 
a few additional measurements made in September 2013 for 
KOI-680. To increase the throughput for these faint targets, the 
high-efficiency mode of SOPHIE was used (resolving power 
R ~ 40 000), as well as the slow read-out mode of the detector. 
The exposure times ranged from 15 to 60 min depending on the 
targets and the weather conditions, in order to keep the signal-to- 
noise ratio as constant as possible for a given target. While the 
first fiber pointed at the target, the sky background was acquired 


1 http://archive.stsci.edu/kepler/data_search/search.php 


Table 1. SOPHIE radial velocity measurements. 


BJD 

-2400000 

RY 

[kms -1 ] 

+ 1 - cr 
[kms -1 ] 

bisect.* 

[kms -1 ] 

exp. 

[s] 

S/Nf 

KOI-206 

55998.6935* 

-1.010 

0.089 

-0.195 

2700 

16.5 

56011.6479 

-1.494 

0.061 

0.062 

2678 

16.4 

56013.6259 

-1.263 

0.062 

0.269 

2685 

13.7 

56037.5652 

-1.486 

0.060 

0.024 

3600 

13.5 

56099.5792 

-1.138 

0.043 

-0.058 

3600 

19.0 

56100.5893 

-1.180 

0.061 

0.351 

2146 

12.4 

56103.4277 

-1.538 

0.043 

-0.252 

3600 

18.9 

56111.5394* 

-1.387 

0.082 

0.029 

3600 

20.6 

56125.5699 

-1.420 

0.047 

0.126 

3600 

16.3 

56153.5630 

-1.127 

0.071 

0.301 

2703 

14.7 

56158.5576 

-1.285 

0.053 

-0.015 

2599 

17.8 

56159.4445 

-1.395 

0.053 

-0.113 

3600 

18.7 

KOI-614 

55830.4663 

-51.355 

0.065 

-0.008 

900 

7.1 

55975.7222* 

-51.096 

0.093 

-0.392 

779 

5.1 

55977.7202* 

-50.994 

0.083 

-0.189 

1126 

5.9 

56038.6270 

-51.252 

0.041 

-0.023 

1800 

9.4 

56042.6242 

-50.933 

0.037 

-0.076 

1800 

12.5 

56086.5856 

-51.350 

0.034 

0.078 

1800 

12.9 

56157.4901 

-50.955 

0.021 

0.073 

1800 

14.7 

56160.5560 

-50.951 

0.029 

0.097 

1800 

11.4 

56161.4765 

-51.088 

0.031 

0.081 

1800 

10.7 

56162.5625 

-51.206 

0.042 

-0.013 

1800 

11.8 

56163.5253 

-51.313 

0.069 

0.074 

1800 

6.7 

KOI-680 

55765.5358* 

-26.361 

0.0238 

-0.0086 

3600 

20.99 

55777.5181 

-26.436 

0.0241 

-0.0087 

2400 

20.62 

55801.5416 

-26.392 

0.0257 

0.0232 

3600 

19.54 

55802.4406 

-26.424 

0.0173 

0.0673 

3600 

29.18 

55804.4901 

-26.451 

0.0205 

0.0398 

2549 

22.90 

55809.4062 

-26.330 

0.0305 

0.1083 

2587 

16.79 

55831.3506 

-26.426 

0.0170 

0.0673 

3600 

28.79 

55832.3516 

-26.347 

0.0182 

0.0489 

3600 

26.47 

55833.3257 

-26.312 

0.0215 

0.0036 

3600 

24.20 

56534.4898 

-26.358 

0.0177 

-0.0426 

3600 

26.90 

56538.4772 

-26.257 

0.0176 

-0.0216 

3600 

29.21 

56552.4366* 

-26.380 

0.0313 

-0.0718 

3600 

20.63 


*: bisector spans; associated error bars are twice those of RVs. 
t: signal-to-noise ratio per pixel at 550 nm. 

J: measurements corrected for moonlight pollution. 


using the second fiber aperture of SOPHIE pointed at ~2 arcmin 
from the target. 


The spectra were extracted using the SOPHIE pipeline 
( [Bouchy et al.|2009] ), which includes localization of the spectral 
orders on the 2D-images, optimal order extraction, cosmic-ray 
rejection, wavelength calibration, and corrections of flat field. 
They were corrected for the spectrograph drifts, measured from 
thorium calibration exposures every two to three hours during 
the night. Radial velocities were obtained from a Gaussian fit to 
the weighted cross-correlation function (CCF) with FO-type nu¬ 
merical masks for the three targets (Baranne e t al.|l996] |Pepe 
et al.||2002} and are reported in Table |jj Some spectra were 
contaminated by moonlig ht and wer e corrected followin g the 
metho d described in Polla cco et al.| ( |2008| ) and |Hebrard et al. 
( 2008) using the background sky obtained with the second fiber. 
The corrected exposures are identified in Table [T] Radial veloc¬ 
ities were also corrected using the RV reference star HD 185144 


( Howard et al.||2010[ |Bouchy et al.||2013[ |Santeme et al.||2014| ) 
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Fig. 1. From left to right: periodograms of KOI-206, KOI-614, and KOI-680. Top panels: Lomb-Scargle periodogram of the com¬ 
plete light curve. Middle panels: mean Lomb-Scargle periodogram of the light curve divided in chunks. Bottom panel: light curve 
autocorrelation. The green dashed lines indicate the center of Gaussians fitted to the data to derive the main frequency and its har¬ 
monics. The red vertical line corresponds to the maximum rotation period deduced from the v sin z* and the star radius, and the 
shaded areas are the associated 1, 2, and 3-cr upper limit regions from the darkest to the lighter ones, respectively. 
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Fig. 2. Bisector span as a function of the detrended radial velocities with 1-cr error bars for KOI-206, KOI-614, and KOI-680 (from 
left to right). Points affected and corrected by Moon contamination are indicated in blue. The ranges here have the same extents in 
the x- and y-axes. 


which was systematically observed the same nights of the three 
Kepler targets. 

The RV dispersion over the two years of measurements on 
this constant star is about 10 m/s but with some variations at the 
level of 20-25 m/s over a few days that seem to be correlated 
with the dome temperature. The correction by this constant star 
slightly improve the residual of the RVs, but without any signifi¬ 
cant change of the parameters. The error bars on the radial veloc¬ 
ities were computed from the CCF using the method presented 
by [Boi sse et al. ( |2010| . For faint targets such as these, some 
spectral domains are too noisy and degrade the accuracy of the 
radial-velocity measurement, therefore the number of spectral 
orders used in the CCF was adjusted to decrease the dispersion 
of the measurements. 

The radial velocity measurements are displayed in the lower 
plots of Fig. [3j together with their best Keplerian fits and 
the residuals. They show variations in phase with the transit 
ephemeris derived from Kepler and with semi-amplitudes im¬ 
plying a companion in the planet-mass range. 

For KOI-614b, no significant linear drift, mask effect, or cor¬ 
relation between the bis sector span (Fig. [2| middle panel) were 
detected, reinforcing the conclusion that the radial-velocity vari¬ 


ations are not caused by spectral-line profile changes attributable 
to blends or stellar activity. On the other hand, the radial veloci¬ 
ties of KOI-206b and KOI-680b show a significant linear drift of 
-500 ± 200 and 43.9 ± 9.2 ms -1 yr -1 , suggesting the presence 
of a distant companion in those two systems of minimum mass 
of about 4.3 ±1.6 and 2.9 ± 0.6 M Jup respectively. 

The nature of these companions is difficult to assess with the 
available data. Using different stellar masks to compute the CCF 
produced no significant effect in both cases. After removing for 
each target the two measurements affected by the moonlight, the 
hypothesis of a null Spearman’s rank correlation coefficient be¬ 
tween the bisector span and radial velocities cannot be rejected 
at the 5% level for both targets, supporting the hypothesis that 
the companion has no detectable effect on the CCF. We find the 
same result when the long-term drift is subtracted from the radial 
velocities as can be seen in Fig. [2] (left and right panels), where 
the CCF bisector spans show neither variations nor trends as a 
function of radial velocity. Neither is there a significative corre¬ 
lation between the radial-velocity long-term drift and the CCF 
bisector spans (Spearman’s rank correlation p-value of 0.49 for 
KOI-206 and 0.32 for KOI-680). 
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The flux of identified (separations > 2") contaminant stars 
in the Kepler photometric mask is estimated for each quarter, 
and the light curves are corrected. There is still the possibility 
that the light curve and the CCF are contaminated by stars closer 
than ~ 2". This contaminant star can be physically related to the 
target or by chance close to the line of sight. We estimated the 
probability of the latter being 1.9, 1.8, and 2.0% for KOI-206, 
KOI-614, and KOI-680, respectively, using the Besangon model 
(Robin et al.] ||2003| ). The magnitude of the stars in the simula¬ 
tion is limited according to the observed transit depth, and its 
distance to the target to less than 2". We try to constrain the 
contaminant adding synthetic spectra to the co-added SOPHIE 
spectrum. The synthetic secondary CCF peak is detectable for 
early than K0V-K3V stars at the same distance of the target (or 
equivalent flux ratio for an unrelated contaminating star), if this 
CCF peak is well detached from the main one in radial velocity. 
That translated into an upper limit for the possible flux contami¬ 
nation in the light curve and CCF of ~ 5% for KOI-206, ~ 10% 
for KOI-614, and ~ 3% for KOI-680. 

Assuming the RV drift observed for KOI-206 and KOI-680 
is caused by a companion in a circular orbit, we estimated the 
allowed masses at each period fitting the RV residuals after hav¬ 
ing subtracted the Keplerian (Section [4]). The result is presented 
in Fig. 11 We estimated the probability that the companion al¬ 
lowed by the drift is stellar in 34% for KOI-206 and 31% for 
KOI-680 using the statistics of Ragh avan et al.| ( [2010| ). As a con¬ 
clusion, for the three targets, the radial-velocity variations de¬ 
tected at the Kepler period are compatible with a companion of 
a planetary nature, and in the case of for KOI-206b and KOI- 
680b, there is no evidence that their estimated parameters are 
significantly affected by the possible additional companions. 


3. Spectroscopic parameters 

For KOI-206 and KOI -680 we used the SOPHIE spectra de¬ 
scribed in Section |Z2| For KOI-614 we used spectra obtained 
with the ESPaDOnS spectrograph. The target was observed in 
‘object+sky’ mode, with a spectral resolution T/AT ~ 65000 
(Program 12AF95). The six observations ran from June 30 to 
July 5, 2012, with exposure times ranging from 3515 s for the 
first one to 1715 s for the five others. Each spectrum was cor¬ 
rected for its radial velocity shift, and then the respective spectra 
were co-added. 

The final spectrum has a signal-to-noise ratio per element 
of resolution at 5500A of 50, 90, and 250 for KOI-614, KOI- 
206, and KOI-680, respectively. We derived the effective tem¬ 
perature T e ff, surface gravity logg, metallicity [Fe/H], and mi 
croturbulence v m i cro of the stars with the VWA package (|Bruntt 
|eFaT. 2004, 2008 


(Gustafsson et al. 


2010b| . We used a grid of MARCS models 
2008| ) to reach the ionization equilibrium over 


a large selection of spectral lines. Then, we fitted a rotational 
profile on a set of isolated metallic lines, solving for the pro¬ 
jected rotational velocity v sin and the macroturbulence v macro . 
For KOI-206 and KOI-680, instead, the v sin /* was measured 
through the Gaussian width of the cross-correlation function 
(CCF) ( [Boisse et al.||2010l and references therein), because of 
the paucity of isolated lines to be fitted with VWA. 

The spectrum of KOI-614 did not present any particular prob¬ 
lem, and the relative atmospheric parameters are presented in 
Table [2] For KOI-206, we found a considerably high value for 
Vmicro (2.8 ± 0.8 km s -1 ). The large uncertainty can be influenced 
by its degeneracy with v macro , which is problematic for con¬ 
straining a spectrum with a low signal-to-noise ratio. However, 


the measurement favors a value higher than 2kms _1 , indicat¬ 
ing an evolved object (e.g., Gray[[2005] ). We imposed different 
values for v macro (between 5.0 and 7.5 kms -1 ) and fitted all 
the other atmospheric parameters. In this way, we verified that 
the latter are not significantly affected. Eventually, we measured 
v sin /*= 11 + 1 km s -1 with the CCF. In conclusion, the ioniza¬ 
tion balance required T& = 6340 ± 140 K, log g = 4.0 ± 0.3 dex, 
and [Fe/H]= 0.06 ±0.19 dex. 

For KOI-680, the analysis yielded T eff = 6090 + 110 K, log g = 
3.5 + 0.1 dex, [Fe/H] = -0.17 + 0.10 dex, and Vmicro = 1.3 + 
0.2 km s _1 . Because of the loose constraint we obtained on v macro 
with VWA, we chose to fix it through the T ef f, according to the 


relation in 


Bruntt et al. 


(2010a), to 3.6kms 

l 


, Finally, from the 


CCF we measured v sin /*= 6 + 1 km s 
We used the SME pack age from | Valenti & Piskunov (1996) and 


Valenti & Fischer ( |2005| ) to check the reliability of the results for 


KOI-206 and KOI-680. We verified T c s on H a , the logg on the 
pres sure-sensitive lines CaT6122, CaT6162, and CaT6439, and 
on the Mglb triplet. Fixing the T e ff, logg, an d [Fe/H] on KOI- 
206, we fitted v macro - 7 km s -1 and v m i cro - 2.8 km s _1 . 

4. Modeling with PASTIS 

We use the Planet Anal ysis and Small Transit Investigation 
Software (PASTIS, Diaz et a l.|2014| ), that allows the light curve, 
the radial velocity, and the spectral energy distribution (SED) to 
be fit with a Markov Chain Monte Carlo (MCMC) algorithm. 
The photometric magnitudes used for the SED fit are listed in 
Table [2] Only the sections of the light curve around the tran¬ 
sits are used, after normalization to a linear trend fitted outside 
the transit. A 3-<x clipping is performed to eliminate outliers, 
and the measured flux is corrected for the contamination factors 
provided for each quarter at MAST archive. PASTIS uses the 
EBOP code ( [Nelson & Davis|1972[|Etzel|1981[|Popper & Etzel 
|1981| ) extracted from the JKTEBOP package ( [Southworth|201 1| ) 
to model the photometric transits. Both long- and short-cadence 
(when available) are used, with an oversampling factor of 10 for 
the long-cadence data to account for the long integration time 
when comparing with the model. Radial velocity curves were 
simultaneously fitted with eccentric Keplerian orbit and a lin¬ 
ear drift. We also account for additional sources of Gaussian 
noise in the light curves, radial velocities, and SED by fitting 
a jitter value to each data set. We used the stellar evolution¬ 
ary tracks from STAREVOL (Palacios, priv. com.), Dartmouth 
( [Potter et al.|2008| ), and PARSEC ( [Bressan et al.|2012~a| to es¬ 
timate the stellar para meters, and the PHO ENIX/BT-Settl syn¬ 
thetic spectral library ( [Allard et al.|[2012| to model the SED. 
The interpolated spectrum is scaled to a given distance and cor¬ 
rected from interstellar extinction. The distance d and the color 
excess E (B _ V ) were included as free parameters in our model. 
The model has therefore 20 jump parameters for KOI-614b, and 
two more for KOI-206b and KOI-680b that have short-cadence 
data, indicated in Table [2] We used non-informative priors (uni¬ 
form or Jeffreys distributions) except for T^, [Fe/H], andp*, for 
which we used spectroscopic estimations as normal and asym¬ 
metric normal distribution priors, and P and T c , for which we 
used the determination by Borucki et al. (2011| as normal distri¬ 


bution priors, but the width of the distribution was increased by 
an order of magnitude to ensure that it would not bias the results. 

It is a well known problem that the uncertainty in the stel¬ 
lar parameters obtained from evolutionary models may be un¬ 
derestimated. To account for this, we use different evolutionary 
models as input for the stellar parameters, and the combined 
outcome will naturally enlarge the errors in the stellar param- 
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Table 2. KOI-206, KOI-614, and KOI-680 system parameters. 


Object 

Kepler-433 

Kepler-434 

Kepler-435 

Kepler object of interest 

KOI-206 

KOI-614 

KOI-680 

Kepler input catalog 

KIC 5728139 

KIC 7368664 

KIC 7529266 

USNO-A2 ID 

1309-0358523 

1329-0421364 

1331-0387694 

2MASS ID 

19502247+4058381 

19342073+4255440 

19290895+4311502 

Coordinates 




RA (J2000) [hh:mm:ss.sss] 

19:50:22.476 

19:34:20.729 

19:29:08.959 

Dec (J2000) [dd:mm:ss.ss] 

Magnitudes 

+40:58:38.17 

+42:55:44.08 

+43:11:50.21 

Kepmag § (AB) 

14.463 

14.517 

13.643 

d51mag § (AB) 

14.745 

14.766 

13.830 

g’ Sloan-Gunn (AB) 

14.953 

14.970 

14.013 

r’ Sloan-Gunn (AB) 

14.413 

14.445 

13.610 

i’ Sloan-Gunn (AB) 

14.253 

14.322 

13.485 

z’ Sloan-Gunn (AB) 

14.166 

14.277 

13.448 

J 2MASS (Vega) 

13.201 + 0.020 

13.400 + 0.022 

12.593 + 0.020 

H 2MASS (Vega) 

12.948 + 0.024 

13.035 + 0.031 

12.345 + 0.019 

Ks 2MASS (Vega) 

12.826 + 0.026 

13.039 + 0.030 

12.291 + 0.022 

WISE W1 (Vega) 

12.716 + 0.025 


12.261 + 0.023 

WISE W2 (Vega) 

12.803 + 0.032 


12.299 + 0.024 

WISE W3 (Vega) 

Spectroscopic parameters 

12.961 + 0.531 



Effective temperature, r eff [K] 

6340 + 140 

5970 + 100 

6090+ 110 

Surface gravity, log g [cgs] 

4.0 + 0.30 

4.22 + 0.10 

3.5 + 0.1 

Metallicity, [Fe/H] [dex] 

0.06 + 0.19 

0.35 + 0.15 

-0.17 + 0.10 

Stellar rotational velocity, v sin /* [ km s -1 ] 

11 + 1 

3.0 + 1.0 

6 + 1 

Microturbulent velocity, v micro [ km s -1 ] 

2.8 + 0.8 

1.0+ 0.2 

1.3 + 0.2 

Macroturbulent velocity, v macro [kms -1 ] 

5 + 1 

2.5 + 1.0 

3.6 (fixed) 

Result from light curve, radial velocity, and SED combined analysis 



Planet orbital period, P [days]* 

5.33408384 + 1.1 x 10 -6 

12.8747099 ± 5.0 x 10“ 6 

8.6001536 ± 1.8 x KT 6 

Mid-transit time, T c [BJD]* 

2454964.98152 ± 1.5 x 1(T 4 

2455003.02283 ± 4.2 x 10" 4 

2455010.64241 ± 4.0 x 10" 4 

cov(P, T c ) [days 2 ] 

-1.25 x 10- 10 

-1.07 x 10" 9 

-2.20 x lO" 10 

Orbital eccentricity, e* 

0.119 ±0.079, <0.29+ 

0.131 +0.072, <0.32 t 

0.114 + 0.077, <0.30 t 

Argument of periastron, oj [deg]* 

68+ m 

89.21+g 

82 + 33 

104 + 36 

Orbit inclination, i [deg]* 

86.46 ^q 38 

85.51 + 0.52 

Orbital semi-major axis, a [AU] 

0.0679 + 0.0027 

0.1143 + 0.0030 

0.0948 + 0.0018 

semi-major axis / radius of the star, a/R+ 

6.44 + 0.62 

17.9 + 1.6 

6.35 + 0.51 

Radius ratio, k = R p /R+* 

0.06590 + 0.00015 

n 

0.06384 + 0.00020 

Linear limb darkening coefficient, u a * 

0.325 + 0.022 

0.00 + 0.43 

0.374 + 0.024 

Quadratic limb darkening coefficient, u b * 

0.220 + 0.045 

0.01 + 0.49 

0.180 + 0.042 

Transit duration [h] 

6.178 + 0.015 

2.288 + 0.045 

9 003 +0 059 
_ 0036 

Impact parameter, b 

0.082 + 0.081 

0.979 + 0.031 

0.448 + 0.024 

Systemic velocity, V r [kms -1 ]* * 

-1.284 + 0.052 

-51.156 + 0.028 

-26.3665 + 0.0079 

Radial velocity linear drift [ms -1 yr -1 ]* 

-500 + 200 

83 + 69 

43.9 + 9.2 

Radial velocity semi-amplitude, K [ms -1 ]* 

259 + 44 

221 + 23 

64+ 12 

Jitter Kepler long cadence, cr Kep i epLC [ppm] * 

65 + 16 

36 + 28 

27 5 +9 - 3 

Jitter Kepler short cadence, cr Kep i er SC [ppm] * 

88 + 58 


48 + 32 

Jitter radial velocity, cr RV [ms -1 ]* 

34 + 32 

16 +19 

0 017 +1 o ° 25 

U ' U1 ' - 0.014 

10 + 10 

Jitter SED, ct S ed [mag]* 

0.056 + 0.024 

0.019 + 0.016 

Effective temperature, T eff [K]* 

6360 + 140 

5977 + 95 

6161 +94 

Metallicity, [Fe/H] [dex]* 

-0.01 + 0.20 

0.25 + 0.14 

-0.18 + 0.11 

Stellar density, p* [p©]* 

0.126 + 0.033 

0.46 + 0.11 

0.047 + 0.011 

Star mass, M* [M 0 ] 

1.46 + 0.17 

1.198 + 0.093 

1.538 + 0.088 

Star radius, R+ [P Q ] 

2.26 + 0.25 

1.38 + 0.13 

3.21 + 0.30 

Deduced stellar surface gravity, log g [cgs] 

3.892 + 0.056 

4.242 + 0.055 

3 613 +0 - 047 

J,u J - 0.070 

Age of the star [Gyr] 

2.67 + 0.91 

4.0 + 1.7 

2.25 + 0.42 

Luminosity of the star, log(L/L 0 ) [dex] 

0.88 + 0.11 

0.33 + 0.11 

1.129 + 0.089 

Distance of the system [pc]* 

1870 + 210 

1240 + 120 

2070 + 200 

Color excess, E^ B _ V ^ [mag]* 

0.210 + 0.039 

0.049 + 0.029 

0.044 + 0.025 

Planet mass, M p [M Jup ] 

2.82 + 0.52 

2.86 + 0.35 

0.84 + 0.15 

Planet radius, P p [R Jup ] 

1.45 + 0.16 

i i O+ 0.26 

1 ’ 1 J - 0.18 

1.99 + 0.18 

Planet mean density, p p [g cm~ 3 ] 

1.13 + 0.32 

2.4 + 1.3 

0.131+0.037 

Planet surface gravity, log g p [cgs] 

3 . 518 !^ 

a 74+O.I3 

7 - 0.21 

2.72 + 0.11 

Planet equilibrium temperature*, T eq [K] 

1776 + 87 

1000 + 45 

1729 + 70 


§ Not used in the SED fitting.* MCMC jump parameter. f upper limit, 99% confidence. * reference time BJD 2456000. 

* T eq = r eff (1 - A) 1/4 using an albedo A = 0. 5 

M q = 1.98842 x 10 30 kg, R Q = 6.95508 x 10 8 m, M Jup = 1.89852 x 10 27 kg, R Jup = 7.1492 x 10 7 m 
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eters from the differences between them. For each of the targets, 
we therefore fit the data three times using a different evolution¬ 
ary model: STAREVOL , Dartmouth, and PARSEC. To ensure a 
broad exploration of the parameter space, we compute 30 chains 
of 10 6 steps, starting at random points drawn from the joint prior. 
After removing the burn-in interval of each chain, we obtain a 
merged chain by thinning each chain. The thinning factor is de¬ 
termined by the maximum correlation length found among all 
parameters in each chain (e.g. [Tegmark et al.|2004] ). This is done 
to obtain independent points to form the posterior distributions. 
Then we take the same number of samples of each merged chain 
from each stellar model, obtaining the combined merged chain. 
We obtain 3168, 729, and 4564 independent points in the com¬ 
bined merged chain for KOI-206, KOI-614, and KOI-680, re¬ 
spectively, which we use to obtain the estimated value and the 
68.3% central confidence interval, both for jump and derived pa¬ 
rameters (listed in Table [2}. The model of maximum likelihood 
is plotted together with the data in Fig. [3] The correlation distri¬ 
butions and histograms are shown in Figs. [12 [13 and 14 for the 
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obtained are solely based on observations. Through the use of 
SET’s MCMC algorithm, we thus obtained the posterior prob¬ 
ability distributions of the bulk composition of the planets (i.e., 
their core mass), as well as an independent estimation of plane¬ 
tary and stellar parameters. The latter are entirely consistent with 
those presented in Section]?] 

The results for KOI-206b, KOI-614b, and KOI-680b are pre¬ 
sented in terms of planetary radii as a function of age in Fig. [6] 
the regions show the 1, 2, and 3-<x constraints from modeling 
the star and transit, while the lines show a subset of planetary 
models for the nominal mass and equilibrium temperature of the 
planet and different compositions. 

All our planetary evolution models assume the following: the 
planet is made of a central rock)0core and a solar-compositioij^] 
envelope. Of course we do not know whether all heavy elements 
are in the core or mixed in the envelope or know their exact com¬ 
position. However, |Baraffe et al. ( [2008j ) show that having the 


MCMC jump parameters. We found that the distributions of the 
posteriors obtained from the different stellar models are compat¬ 
ible within 1-cr. The stellar evolution tracks from STAREVOL 
in a Hertzsprung-Russell diagram are given in Fig. [4] along with 
the luminosity-r ef f joint posterior distributions to show the evo¬ 
lutionary stage of the host stars. 


5. Planetary evolution models 

KOI-614b, KOI-206b, and KOI-680b are all different from each 
other: with radii of 1.13, 1.45, and up to 1.99 Rj up , and masses of 
2.86, 2.82, and down to 0.84 M Jup , they occupy distinct regions 
of the mass-radius diagram as illustrated in Fig. [5] 

To estimate the structure and composition of the three plan¬ 
ets, we computed the joint evolution of b oth planet and host star 
for each syste m, combining t he CEPAM ([Guillot & Morel|l 995 [ 
Guillot|20"l0] ) and PA RSEC ([Bressan et al.|2012b|) evolutionary 
models using SET (see|Guillot & Havel|2011[|Havel et al.|2011| 
see also Almenar a et al.|2013| ). Since stellar and planetary abso¬ 
lute parameters are ultimately model dependent, the values we 


heavy elements mixed in the envelope generally leads to smaller 
radii. Therefore our models should provide an upper limit for Z, 
the fraction of heavy elements in the planet. In addition, we do 
not expect the envelope to have a solar composition (it is not even 
the case for Jupiter or Saturn; see, e.g., [Guillot & Gautier|2014| ), 
but the uncertainties in the models themselves (EOS, opacities, 
atmospheric model) play a larger role than the small adjust¬ 
ment^ in the composition of the envelope due to the different 
stellar metallicities (Z* ~ 0.013 for KOI-680 and Z* ~ 0.042 
for KOI-614). 

For consistency with previous modeling of exoplanets, plan¬ 
etary evolution models are calculated in two cases: (i) a “stan¬ 
dard” case where only the irradiation from the star is accounted 
for; (ii) a “dissipated-energy” model in which, in addition to the 
standard case, a fraction (0.25%) of the incoming stellar flux is 
assumed to be converted into kinetic energ)p]and thendissipated 
at the center of the planet (see, e.g., |Guiliot & Showman 2002 


and Spiegel & Burro wsf20 13 1 for detailed discussions). 

With a radius of 1.13 Rj up and a mass of 2.86 Mj up , both 
standard and dissipated models provide solutions for the plan¬ 
etary radius of KOI-614b that match the available constraints, 
resulting in a core mass of 81^ 4 M®. As seen in Fig.p] models 
with no core (i.e., pure H/He planet) land in the miacrfe of the 
1-cr region, but the high metallicity of the star (+0.35) tends to 
indicate that the planet should probably h ave a significant num¬ 
ber of heavy elements ( [Guillot et al.|2006] ). 

On the other hand, KOI-206b has a similar mass of 
2.82 Mj up , but a radius of 1.45 Rj up , and an inferred density of 
1.13 gem -3 . For such a high radius and nearly three times the 
mass of Jupiter, KOI-206b does not have many similar siblings: 
the t hree closest in the m ass-radius diagram (Fig. [5]) are CoRoT- 
2b ( [Alonso et al.|2008|), CoRoT-18b ([Hebrard et al.|2011| ), and 


Kepler-17b (Pese rTet al.|201 1[ Bonomo et al. 2012). The former 
two, for instance, have been modeled and are still challenging 
the current_planetary models (see, e.g., [Guillot & HaveT||2011[ 


Mou tou et al.||2013] ). As Fig. [6] shows, this is also the case of 
KOI-206b: standard models barely succeed at matching the 2- 
cr age-radius constraint, while dissipated-energy models are just 
reaching the bottom of the 1-cr region. Although an indepen- 


Fig.5. Radius versus mass of all known transiting planets as of 
February 2014 (see http://exoplanets.org/, [Wright et al.pOll) ). 
The three planets presented in this paper are highlighted in red 
and labeled. 


2 38% Si0 2 , 25% MgO, 25% FeS, 12% FeO 

3 70% H, 28% He, 2% others (Z) 

4 as c omp ared t o the exp ected, much higher Z values in planets than 
in stars ([Guillot et al.|2006|) 

5 Note the use of ‘ke’ labels in figures: 1% ke = 0.25% of the incident 
flux at substellar point or 1% of the irradiating flux averaged over the 
planet’s area converted into kinetic energy (ke). 
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KOI-614 


KOI-680 


Wavelength [A] 
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BJD - 2,400,000 


56400 56600 




Fig. 3. Model of the maximum likelihood step of the merged chain of KOI-206, KOI-614, and KOI-680 (from left to right columns) 
plotted with the data. From top to bottom: SED, light curve, radial velocity in time, and radial velocity in phase. In the bottom 
panels: residuals after subtracting the model to the observed data. SED plot: the solid line plots the PHOENIX/BT-Settl interpolated 
synthetic spectrum, in red circles the absolute photometric observations listed in Table [2J in blue open circles the result of integrated 
the synthetic spectrum in the observed bandpasses. Light curve plot: folded transit observed by Kepler in long-the black solid line 
plots the Keplerian model a for the radial velocities and red circles the SOPHIE observations. 





T eS [K] T eff [K] T eff [K] 

Fig. 4. Stellar evolution tracks from STAREVOL , from the zero-age main-sequence stage to 25 Gyr age, in the luminosity-effective 
temperature plane for [Fe/H] = 0.00, 0.25, and -0.20 (from left to right, solid line) and [Fe/H] = 0.10, 0.35, and -0.15 (from left to 
right, dashed line) with the joint posterior distributions of the MCMC (98.9% joint confidence region is denoted in red) for KOI-206, 
KOI-614, and KOI-680, respectively. The mass in solar masses is annotated at the beginning of the main sequence of each track. 
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0 2 4 6 8 10 12 14 

age [Ga] 


Fig. 6. Evolution of the planet radius as a function of the age. KOI-206b with mass 2.82 Mj up and effective temperature 1776 K is 
plotted to the left, KOI-614b (2.86 M Jup and 1000 K) in the center plot, and KOI-680b (0.84 M Jup and 1729 K) in the right plot. The 
68.3%, 95.5%, and 99.7% joint confidence regions are denoted by black, dark gray, and light gray areas, respectively. The curves 
represent the thermal evolution of a respective planet mass and equilibrium temperature. Text labels indicate the amount of heavy 
elements in the planet (its core mass in Earth masses). Dashed lines represent planetary evolution models for which 0.25% of the 
incoming stellar flux (e) is dissipated into the core of the planet, whereas plain lines do not account for this dissipation (standard 
models). Dashed-dotted lines (gray shades) are coreless models with higher amounts of dissipated energy, as labeled. 


dent MCMC 1-D distribution gives a core mass of 0 M® with 
a 68.3% confidence interval [0, 74] M®, thus promoting solu¬ 
tions for which the planet has a much lower radius. Extra mod¬ 
els with incredibly large amounts of energy dissipated in the 
planet have to be created to explain the observed radius at 1- 
<T. Up to 75% of the total light reaching the planet would need 
to be dissipated for these models to work out a solution. There 
is, to date, no known physical proces s that would allow such ef- 
ficiency. However, the recent work of Spiegel & Burrows ( 2013) 
suggests that dissipating the energy in the upper atmosphere in¬ 
stead of the deep interior maintains the planet’s “hot state” (i.e., 
large radius) longer for ages below 2-3 Gyr. Last but not least, 
the use of a newer H/He equation-of-state (EOS) may contribute 
significantly to solving the problem: while our models are using 
SCVH (S aumon-Chabrier-Van Horn, |Saumon et al.||1995] ), the 
EOS from Militzer & Hubbard ( 2013]) may produce planets with 


bigger radius, up to about 0.2 R Jup larger than ours. 

Last of all, with a radius of 1.99 Rj up and a mass of 0.84 M Jup 
KOI-680b is among the largest planets ever kn own. Onl y tw o 
planets have a higher radius, HAT-P-32b ([Hartman et a l.|2011| ) 
and WASP-17 ( Anderson et al. 201 1| ) (see Fig. [5]), and as few as 


eleven have a radius larger than 1.7 Rj up . WASP-78b (Smalley 
et al.|2012| ) is the only one of similar mass and radius around an 


evolved star with low metallicity: unfortunately, it has not been 
modeled yet. As expected, our standard planetary models cannot 
account for such a high radius. Even with the dissipation of 1% 
of the available irradiation, we hardly reach the 3-cr constraint 
region. Once again, we had to produce some more extreme mod¬ 
els with unphysical amounts of energy dissipated into the deep 
planet to achieve the 2 R Jup limit. Interestingly, we note that dis¬ 
sipating “just” 10% to 40% of the incident stellar flux is enough 
to match the 1-cr region, possibly less than what is required for 
KOI-206b. Previously suggested improvements over our models 
for KOI-206b should apply for this planet. On the other hand, 
KOI-680b is below a Jupiter mass and should therefore be easier 
to re-inflate. 


6. Discussion 

6.1. KOI-614b 

KOI-614b was a moderate probability candidate, but the ra¬ 
dial velocity measurements presented in this paper establish the 



Fig. 7. Planetary mass as a function of the orbital period for 
known transiting giant exoplanets |Wright et al.|201 1| ). The three 
planets presented in this paper are highlighted in red and labeled. 
The dotted line represents the Roche limit, and the dashed line 
corresponds to twice this value, which would be the final orbital 
period of an initial eccentric long-period planet that would have 
completed its orbit. 


planetary nature of the companion, its high-impact parameter 
(b « 0.98) explaining the V-shaped transits. With an equilibrium 
temperature of T eq = 1000 ± 45 K, KOI-614b is at the transi¬ 
tion between what is usually referred to as “hot” and “warm” 
Jupiters and is a good candidate for testing migration theories. 
KOI-614b is one of the few transiting giant planets with an or¬ 
bital period longer than ten days, as can be seen in Fig. [ 7 ] Unlike 
hot Jupiters, this kind of planets has a pericenter distance that is 
too large to allow efficient tidal dissipation to induce migration 
through the high-eccentricity migration mechanism ( |Dong et al. 


2014). In order to access the close pericenter required for mi¬ 


gration and circularization by tidal dissipation of the eccentric¬ 
ity, they need to be accompanied by a strong enough perturber 
to overcome the precession caused by General Relativity. KOI- 
614b has a low eccentricity (e < 0.32), but there is no trace of 
an additional massive companion in the system. Using the radial 
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greater heavy element content can be obtained with particular 
disk conditions (the case of “monarchical growth”, Mordasini 
[et al.||2009t |Weidenschilling||2005| ) or by the collision with an 


Fig. 8. Planetary density as a function of mass for known tran¬ 
siting giant exoplanets ( [Wright et al.|||2011 ). The three planets 
presented in this paper are highlighted in red and labeled. 


velocity residuals of our fit, we can exclude the presence of a 
companion of more than 2.0 M Jup on orbital periods less than 
200 days, at the 3-cr level. Given the length of our observation 
period, we cannot make any conclusion about more distant com¬ 
panions. In addition, Mazeh et al. ( |2013| ) report no significant 
transit timing variations for this target. The high-eccentricity mi¬ 
gration scenario cannot be favored yet. 

Considering the mass of - 3 Mj up of KOI-614b, its for¬ 
mation would have quickly opened a gap in the protoplanetary 
disk, and another migration scenario would involve type II mi¬ 
gration. However, as shown by [Burkert & Ida[(2007| ), whatever 
the disk lifetime, type II migration of a single Jovian planet in 
a truncated disk naturally leads to the creation of period gap be¬ 
tween 0.08 < a < 0.6 AU in the distribution semi-major axis, 
which is especially pronounced for systems with stellar masses 
M* > 1.2 M 0 . As this is precisely the case of this system, 
the migration of the planet within the disk could be doubted. 
Nevertheless, while single planets open only narrow gaps, mul¬ 
tiple planets can clear wide inner holes ( [Dodson-Robinson ~&] 
Salyk|20lT] ), thus preventing the migration of the Jovian planet 


down to a short-period orbit. As previously said, our observa¬ 
tions exclude the presence of Jovian companions over orbital 
periods <100 days, but we cannot exclude the presence of less 
massive, terrestrial planets in the system. In this regard, the com¬ 
position of the planet can provide additional clues. The estimated 
radius of - 1.13 Rj up makes KOI-614b a high-density giant 
planet. As can be seen in Fig [8] this mass and density are com¬ 
mon in the population of the transiting exoplanets discovered 
so fai[] Within our estimated uncertainties, we have established 
that the dimensions of this object can encompass internal com¬ 
positions ranging from pure hydrogen to containing a solid core 
of heavy elements up to - 150 M®. Even if the pure hydrogen 
scenario cannot be discarded, some authors (see, e.g., |Guillot| 
|et al.|2006l > have noted that there is observational evidence of a 
correlation between star metallicity and heavy elements in the 
planet. Considering that KOI-614b is orbiting a main-sequence, 
metallic ([Fe/H] = 0.35 dex), GOV star, it is reasonable to fa¬ 
vor a high heavy-element content. In the core-accretion scenario, 
the “canonical” core mass is about 10 M® ( |Mizuno|[l980| . A 


as of Februray 2014 at http://www.exoplanets.org 


other planet(s) or planetesimals after substantial accretion of the 
envelope (Ikoma et al. 2006[ ). In the latter case, the disk would 
have fostered the formation of multiple embryos, which would 
provide a natural way to stop the migration of the planet at its 
current semi-major axis («a -0.114 AU). 

Both migration scenarios require the presence of additional 
companion(s) in the system. High-eccentricity migration mech¬ 
anisms require a massive, not very distant perturber to bring the 
warm Jupiter into a Kozai-Lidov cycle, type II migration re¬ 
quires multiple planets to clear a wide inner hole and prema¬ 
turely stop the migration of the planet. KOI-614b is thus an in¬ 
teresting candidate for future observations with the aim of find¬ 
ing additional low-mass, short-period planetary companions or 
a more distant and massive perturber. 

6.2. KOI-206 

Like KOI-614b, KOI-206b is a relatively massive planet of about 
2.8 M Jup , but its radius of 1.45 + 0.16 Rj up implies a different 
inner structure with a Jupiter-like density. Although the much 
tighter orbit of KOI-206b (<a - 0.0679 AU) puts it in the range 
of hot Jupiters, and although the higher effective temperature of 
its host star does imply a stronger irradiation than for KOI-614b, 
its large radius is particularly challenging for planetary evolution 
models. 

The radius of a transiting planet is one of the most robust 
measurements because it only relies on the knowledge of the 
transit depth and of the stellar radius (with a second-order de¬ 
pendence on limb darkening, and perhaps light curve contami¬ 
nation). Here, the stellar radius is determined using evolution¬ 
ary models, which are known to be inaccurate, and yet this is 
partially taken into account by our fitting method by combining 
several sets of evolutionary tracks in the final result. Concerning 
the transit depth, it is precisely measured at the 210-cr level. 
Nonetheless, the light curve of KQI-206 b ex hibits rotational 
modulations due to spot coverage (cf. Sect. |2.1| ). 

There is no sign of spot occultation during the transits of the 
planet, but high-latitude spot could still affect the estimation of 
the radius. In particular, the presence of such a group of spots 
could be indistinguishable in the light curve and would induce 
an overestimation of the radius of the planet. As shown in Fig.[6j 
typical models of bloated hot Jupiters, where 0.25% of the in¬ 
coming stellar flux is dissipated in the core, would agree better 
with a planetary radius in the bottom of the l-cr confidence re¬ 
gion. Still, to reduce the estimated planetary radius to this value, 
i.e. to about 1.3 Rj up , we would need to assume that the value of 
the flux continuum is actually a factor 1.26 more than the value 
we are using. It would require a polar spot of angular radius of 
25°, in the very conservative assumption that the contrast of the 
spot is 0 and the rotational axis of the star is aligned on the line 
of sight. 

More realisticall y, if we assume a surface-to-spot contrast of 
about 0.67 ( [Lanza et~aL 2003) and the limb darkening coeffi¬ 
cients of our best model, this would req uire a polar spo t with 
a radius of at least 59° (estimated using [Kipping| 2012| ). Such 
a spot could be betrayed by the activity level estimated with 
the emission cores of the Call lines. Unfortunately, the signal- 
to-noise ratio of the SOPHIE spectra in the blue part of the 
spectrum is too low to estimate the activity with Call H and K 
lines. Either way, the dimensions of such a spot is rarely ob¬ 
served ( Lanza et al.|2003 ), and it is very unlikely that our plan- 
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Star radius [R®] 


Fig. 9. Planetary radius as a function of stellar radius for known 
transiting exoplanets ( [Wright et al.|201l| . The dotted lines give 
the loci of some typical transit depth given in labels. The three 
planets presented in this paper are highlighted in red and labeled. 


etary radius estimation would be overestimated because of such 
an effect. Besides, the stellar density inferred from the transit 
(Table. [2} is in good agreement with the stellar parameter deter¬ 
mined by the spectral analysis, and there is no evidence of a false 
positive. 

At the time of writing, KOI-206b is the transiting planet de¬ 
tected around the host with the second largest radius so faiQjust 
after KOI-680b, as shown in Fig. [9] There are only a few other 
known bloated planets in this mass range (see Fig. [5ji, and KOI- 
206b is the one with the lowest density in this group (see Fig.[8]>. 
For such a mass, the amount of extra dissipated energy to pro¬ 
duce that kind of density is preposterous. KOI-206b orbits an 
evolved F7-type star that has recently left the main sequence (see 
Fig.0. and the dynamic aspects related to the rapid and recent 
evolution of the host might have effects that are unaccounted for. 
As noted by |Burrows et al.| ( 2003| ), there is a difference between 
the observed transit radius (at -rnbar levels) and the modeled 
radius, AR. Since the evolution of a giant planet is largely con¬ 
trolled by its atmosphere, the planetary evolution code we used, 
CEPAM, has a dedicated module that computes the exact P-T 



ing the interior structure of the three giant planets can be directly 
compared to the observations: AR cannot be the explanation for 
the discrepancy between theoretical models and observed radius. 


6.3. KOI-680b 

The host star of KOI-680b shares many properties with the one 
of KOI-206b, but with a mass of only 0.84 ±0.15 M Jup and a ra¬ 
dius of 1.99 + 0.18 Rjup, KOI-680b is one of the biggest (Fig.|5]) 
and least dense (Fig. [8]) planets detected to date. Such properties 
exclude a high concentration of heavy elements in its interior 
(Fig.0. which correlates well with the low metallicity of the 
host star ([Fe/H] = -0.17 ± 0.10 dex). Although standard plan¬ 
etary evolution models fail to reproduce its observed character¬ 
istics, the low mass of KOI-680b also means that it is easier to 

7 see note [6] 
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Fig. 10. Kiel diagram for the stars of known transiting (in blue) 
and radial velocity (in orange) exoplanets ( [Wright et al.|2011 ). 
KOI-614, KOI-206, and KOI-680 are drawn in black and la¬ 
beled. We also plot the STAREVOL evolutionary tracks for stars 
of 1, 1.5, and 2 M 0 and solar metallicity. 


inflate, and improv ements in the equation of state might easily 
solve the problem ( jMilitzer & Hubbard||2013| . Not only is the 
planet one of the biggest ever observed, but it also orbits the 
bigger star observed with the transit technique (Fig. [9]). Indeed 
the host star of KOI-680b is an F9 subgiant, well on its way to¬ 
ward the beginning of the red giant branch (RGB). Although it 
is slightly younger than the host of KOI-206b, it is also more 
massive (M* = 1.538 ± 0.088 M 0 ) and thus more evolved. As 
can be seen in Fig. [TO] the host star of KOI-206b is observed in 
a phase of rapid stellar evolution, where no other planets have 
been detected before. 

With an orbital semi-major axis of 0.0948+0.0018 AU, KOI- 
680b is doomed to be engulfed by the star, even neglecting tidal 
interaction, the stellar radius will catch up with the current orbit 
of the planet within the next 260 Myr or so. We note that we de¬ 
tect a companion of minimum mass 2.9+0.6 M Jup with minimum 
orbital period of about 790 days (i.e., about 1.9 AU). Further ra¬ 
dial velocity follow-up is required to determine the period and 
the planetary nature of the companion. 

7. Conclusion 

We have presented the detection and characterization of KOI- 
614b, KOI-206b, and KOI-680b, which are three new transit¬ 
ing, giant extrasolar planets. They were first detected as promis¬ 
ing candidates by the Kepler team from the Kepler light curve. 
We established their planetary nature with radial velocity follow¬ 
up using the spectrograph SOPHIE and characterized the sys¬ 
tems through a combined MCMC Bayesian fit of photometric 
and spectroscopic datasets. We find that KOI-614b is a dense, 
warm Jupiter, with a mass of 2.86 + 0.35 Mj up and a radius of 
l-13^oi8 Rjup, and it is orbiting a GO, metallic ([Fe/H]=0.35 + 
0.15) dwarf in 12.9 days. KOI-206b has a mass of 2.82 + 
0.52 M Jup and a radius of 1.45 + 0.16 Rj up , and it orbits a slightly 
evolved F7-type star in a 5.3-day orbit. Orbiting a F9 subgiant in 
8.6 days, KOI-680b has a much lower mass of 0.84 + 0.5 M Jup 
and an uncommonly large radius of 1.99 + 0.18 R Jup . With host 
stars of masses of 1.46 + 0.17 M 0 and 1.54 + 0.09 M© and 
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radii of 2.26 ± 0.25 R Q and 3.21 ± 0.30 R Q , KOI-206b and KOI- 
680b are the transiting planets around the two biggest stars ob¬ 
served so far. They are inflated hot Jupiters, which are partic¬ 
ularly challenging for planetary models that require unusually 
large amounts of additional dissipated energy in the planet, es¬ 
pecially the massive KOI-206b. For those two targets we also 
find signs of a distant companion that are compatible with the 
planetary mass range. Although we did not perform a detailed 
analysis of the false positive scenarios, we note that there is no 
evidence of a blend or of unaccounted-for contamination. There 
is no significant spatial shift at the time of the transits in the 
Kepler photometric centroid. The transits of KOI-206b and KOI- 
680b are well sampled and U-shaped, with a large number of 
points in the ingress and egress phases. There are no detectable 
lines of a contaminant in the spectra or in the CCF that shows 
any mask effect. The bisector span is not correlated to the radial 
velocities at the 95% level. For both targets, the projected rota¬ 
tional velocity of the star is compatible with the modulations of 
the light curves due to spots. The stellar density determined with 
the transit is in good agreement with the one determined by the 
spectroscopic analysis. Using three different stellar evolutionary 
tracks, the inferred stellar masses, radii, and ages are compatible 
within 1-cr. Nevertheless, more observations of those systems 
would be interesting for assessing the nature and orbital period 
of the companion. That would provide new insight into the sur¬ 
vival of planets through the late evolution stages of their stars 
and help improve the theories of formation and survival of the 
short-period planets around massive stars. 
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Fig. 11. Third-body mass allowed by the RV drift for a given period assuming circular orbits for KOI-206 (left) and KOI-680 (right). 
The median of the distribution is represented by a solid line, and the 1-, 2-, and 3-<x confidence intervals are represented with three 
different gray levels. The dotted lines separate the planetary, brown dwarf, and star regimes. 
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Fig. 12. Two-parameter joint posterior distributions of the MCMC jump parameters of the fit with stellar models of KOI-206. The 
39.3%, 86.5%, and 98.9% joint confidence regions are denoted by three different gray levels. The histogram of each parameter is 
shown at the top of each column, except for the parameter on the last line that is shown at the end of the line. 
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Fig. 13. Idem Fig. 
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but for KOI-614. 
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but for KOI-680. 





































































